Breaking rate minimum predicts the collapse point of over-loaded materials 
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As a model of composite materials, we choose a bundle of fibers with stochastically distributed 
breaking thresholds for the individual fibers. The fibers are assumed to share the load equally, and 
to obey Hookean elasticity right up to the breaking point. We study the evolution of the fiber 
breaking rate at a constant load in excess of the critical load. The analysis shows that the breaking 
rate reaches a minimum when the system is half-way from its complete collapse. 

PACS numbers: 02.50.-r 
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I. INTRODUCTION 



Bundles of fibers, with statistical distributed thresh- 
olds for the breakdown of individual fibers, present inter- 
esting models of failures in materials. They have simple 
geometry and clear-cut rules for how stress caused by 
a failed element is redistributed on undamaged fibers. 
Since these models can be analyzed to an extent that is 
not possible for more complex materials, they have been 
much studied (For reviews, see [3, S H, BL HI ) • The statis- 
tical distribution of the size of avalanches in fiber bundles 
is well studied [1, 0, G3 > an d the failure dynamics under 
constant load has been formulated [|[ through recursion 
relations which in turn explore the phase transitions and 
associated critical behavior in these models. 

In this article we present a way to predict when an 
over-loaded bundle collapses, by monitoring the fiber 
breaking rate. 

We focus on the equal-load-sharing models, in which 
the load previously carried by a failed fiber is shared 
equally by all the remaining intact fibers [TJ E2> EH E3 • 
We consider a bundle consisting of a large number N of 
elastic fibers, clamped at both ends (Fig. 1). The fibers 
obey Hooke's law with force constant set to unity for 
simplicity. Each fiber i is associated with a breakdown 
threshold Xi for its elongation. When the length exceeds 
Xi the fiber breaks immediately, and does not contribute 
to the strength of the bundle thereafter. The individ- 
ual thresholds assumed to be independent random 
variables with the same cumulative distribution function 
P{x) and a corresponding density function p(x): 



Prob(iEi < x) = P(x) 



p(y) dy- 



(1) 
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FIG. 1. The fiber bundle model. 

If an external load F is applied to a fiber bundle, the 
resulting failure events can be seen as a sequential process 
[8|. In the first step all fibers that cannot withstand 
the applied load break. Then the stress is redistributed 
on the surviving fibers, which compels further fibers to 
fail, etc. This iterative process continues until all fibers 
fail, or an equilibrium situation with a nonzero bundle 
strength is reached. Since the number of fibers is finite, 
the number of steps, tf , in this sequential process is finite. 

At a force (or elongation) x per surviving fiber the total 
force on the bundle is x times the number of intact fibers. 
The expected or average force at this stage is therefore 



F(x) = Nx(l- P(x)). 



(2) 



The maximum F c of F(x) corresponds to the value x c for 
which dF/dx vanishes. Thus 



1 - P(x c ) - x c p(x c ) = 0. 



(3) 



We characterize the state of the bundle as pre- critical or 
post-critical depending upon the stress value a = F/N 
relative to the critical stress 
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Fc/N, 



(4) 
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We study the stepwise failure process in the bundle, when 
a fixed external load F = Na is applied. Let N t be the 
number of intact fibers at step no. t, with No = N. We 
want to determine how Nt decreases until the degrada- 
tion process stops. With N t intact fibers, an expected 
number 



[NP(Na/N t )} 



(5) 



of fibers will have thresholds that cannot withstand the 
load, and consequently these fibers break immediately. 
Here [X] denotes the largest integer not exceeding X. 
The number of intact fibers in the next step is therefore 



JVt+i = N — [NP(Na/N t )} 



Since N is a large number, the ratio 



(6) 



n t 



Nt 
N 



(7) 



can for most purposes be considered a continuous vari- 
able. By © we have essentially @ 



n t+ i = 1 - P{a/n t ). 



(8) 



II. THE RELATION BETWEEN MINIMUM 
BREAKING RATE AND COMPLETE COLLAPSE 



We will now demonstrate, for three different threshold 
distributions, that there is a a relation between the min- 
imum of the breaking rate R(t) = —drit/dt (treating t as 
continuous) and the moment t / when the complete fiber 
bundle collapses. 



Uniform distribution 



We consider the uniform distribution, P(x) — x for 
< x < 1, and assume that the load is post-critical: 
a = j + e, with e > 0. Simulations show that the 
breaking rate has a minimum at some value io(e), an d 
that for varying e the minima all occur at a value close 
to | when plotted as function of the scaled variable t/tf 
(Fig. 2). 




FIG. 2. The breaking rate R(t) vs. step t (upper plot) 
and vs. the rescaled step variable tf/t (lower plot) for the 
uniform threshold distribution for a bundle of iV = 10 7 fibers. 
Different symbols are used for different excess stress levels 
a — a c : 0.001 (circles), 0.003 (triangles), 0.005 (squares) and 
0.007 (crosses). 

This can be shown analytically. The iteration ((5|) takes 
in this case the form 



nt+i 



(i + e) 



1 



(9) 



By direct insertion one verifies that 

nt = \ — \fet&n(At — B), (10) 

where 

A = taja" 1 (2 v / e) and B = tan _1 (l/2v / e), 

is the solution ([9]) satisfying the initial condition no = 1. 
From (JTOJ) follows the breaking rate 

R(t) = = yfi A cos~ 2 ( At- B). (12) 

i?(£)has a minimum when 
dR 







dt 



cx sin(2At - 2B), 



(13) 
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which corresponds to 



B 

t = T 



(14) 



When criticality is approached, i.e. when e — ► 0, we have 
A — > 0, and thus t — > oo, as expected. 
We see from eg. (ITU|) that n t = for 



(B + tan" 1 (1/2 Vi) /A = 2B/A. 



(15) 



This is an excellent approximation to the integer value 
at which the fiber bundle collapses completely. 

Thus with very good approximation we have the simple 
connection 



tf = 2tc, 



(16) 



When the breaking rate starts increasing we are halfway 
to complete collapse! 



B. Displaced uniform distribution 



Consider a uniform distribution on the interval (xi, 1): 



■ xi < x < 1 

p(x) = < 1 ~ x t ~ ~ 
' ' otherwise 



Thus 



P(x) 



x < xi 

1— X; 



Xi < X < 1 



(17) 



(18) 



Simulations of the breaking rate gives qualitatively 
the same behavior as for the uniform distribution (Fig. 3). 




FIG. 3. The breaking rate R(t) vs. the rescaled step variable 
tf/t for the displaced uniform threshold distribution (117[1 . 
Here xi = 0.2 and N = 5 x 10 7 . Different symbols are used 
for different excess stress levels a ~ a c : 0.001 (circles), 0.003 



(triangles), 0.005 (squares) and 0.007 (crosses) 
For this distribution eq.@ gives 

x(l — x) 



a = x(l - P(x)) 



1 - xi 



with a maximum 



1 



4(1 -or;) 



(19) 



(20) 



at x — x c — 1/2. 

The iteration (El) takes now the form 



«t+i 



1 - a/n t 



(21) 



1 - xi 

This can be cast in a familiar form. Introduce 

y t = n t (1 - xi) 
in (j2"Tj) to obtain the iteration for y t : 

Vt+i = 1 - c(l - xi) ■ — . 

Vt 

By (|20l) the critical value of cr(l — x;) is 1/4, so we may 
write 



(22) 



(23) 



<j(l - X;) 



(24) 



where again e is assumed to be small and positive. Then 
we are back to the same iteration ^ as for the usual 
uniform distribution: 



y t = h - \fl tan 



tan 



-1 ( 2 



-^-j +t tan- 1 (2Vi) 



(25) 



or, since y = n (l — xi) = 1 - xf. 

yt = g-Ve tan[-tan" 1 ((l/2-Xi)/\/e)+ttan" 1 (2\/i)] 

(26) 

For simplicity write this as 

yt = \ - Ve taa(at - 6), (27) 

with 

a = tan _1 (2 v / e) and 6 = tan _1 ((l/2 - xi)/y/e) 

(28) 

The breaking rate (treating t as continuous) is 



DM dn t _Ldy t yg 



1 — IE; dt 1 — Xi 



cos 2 (af - b) 



(29) 



The minimum breaking rate occurs when dR/dt oc 
sin(2ai — 2b) = 0, i.e. at to = b/a. For small e we use 
the identity 



tan _1 (l/r/) = it/2 - taiT 1 ^), 



(30) 
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and obtain approximately for small e 

a ps and b ps tt/2 - y/e/(l/2 - xi) 
Using this, we obtain to leading order 

to 



7T 



(31) 



(32) 



A good approximation to the collapse point t/ is ob- 
tained by selecting the t for which n t or y t vanishes. From 
(|27|) we see that this occurs for at/ given by 



l~i tan(at/ — 6) = 0, 



= [6 + tan _1 (l/2\/e)]/a. 
Again, by using ([30| we have for small e 
tt/2 - y/i/(l/2 - xi) + tt/2 - 2^1 



(33) 
(34) 



(squares) and 0.007 (crosses). 

This case is more complicated, but the analytical 
ground work has already been done in Q . Eq. (29) in 3 
shows that for small e the iteration is of the form 



(37) 



n t — n c — byJe/C tan(tv Ce — c) 



Here n c = e 1 / 5 , C — l^e) 1 / 5 , b = 5 1 / 5 , and the con- 
stant c is determined by the initial condition hq = 1: 



c = tan 



(l-n c )6^ i yCA 



(38) 



From (f3"T)) . the breaking rate equals 

R(t) = ~^-ot cos- 2 {tVC~e~c). (39) 



Comparing the results for to and t f we have once more 
tf/h = 2. (36) 

to leading order. 



_ _ The breaking rate is a minimum when the cousins takes 

— —= (l + C(\/e))its maximum value 1. This is the case when 
e 



tn = 



(Ce) 



-1/2 



tan 



{l-n^y/C/e . (40) 



The inverse tangent is close to 7r/2 when e is very small. 
Hence, for small overloads, we have in excellent approxi- 
mation 



tn 



C. Weibull distribution 

Let us finally consider a completely different thresh- 
old distribution, a Weibull distribution of index 5, 
P{x) = 1 — e~ x \ Simulations reveal that the breaking 
rate has a similar behavior as in the two cases considered 
above (Fig. 4). 




FIG. 4. The breaking rate R(t) vs. the rescaled step variable 
t/tf for a bundle of N — 10 7 fibers having a Weibull threshold 
distribution. Different symbols are used for different excess 
stress levels a — cr c : 0.001 (circles), 0.003 (triangles), 0.005 



2VCe 



(41) 



The collapse point tf is already evaluated in Q, with 
the result 



(42) 



Cc 



for small e (eq. (33) in 0]). 

Comparison between (|42[) and (|4"Tj) gives 



t 



f 



2t, 



(43) 



as for the two previous threshold distributions consid- 
ered. 



III. COMMENTS 

We have shown that the complete collapse of fiber 
bundles occurs at tf = 2to, where to denotes the number 
of steps of the breaking process at which the fiber 
breaking rate has a minimum. The results are derived 
for very small overloads e. For larger overloads the ratio 
to /tf will not be exactly 0.5, as illustrated in Fig. 5, but 
nevertheless of the order of 0.5. 
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FIG. 5. Simulation results for the ratio to/tf vs. (a — <7 C )~ 1 '' 2 
for the uniform distribution (circles), the displaced uniform 
distribution with Xi = 0.2 (squares) and for the Weibull 
distribution (triangles). The graphs are based on 1000 
samples with N = 10 7 fibers. 



FIG. 6. Simulation results for n(to) vs. (a — a c ) for the uni- 
form distribution (circles), the displaced uniform distribution 
with xi — 0.2 (squares) and for the Weibull distribution 
(triangles). The graphs are based on 1000 samples with 
N = 10 7 fibers. The straight lines represent the critical value 
n c for these three distributions. 



Another interesting observation is that at t = to the 
number of unbroken fibers in the bundle n(trj) attains 
the critical value n c . This can be derived analytically by 
putting the value of to in the expressions (10), (27) and 
(37) respectively, for the uniform, the displaced uniform 
and the Weibull distribution. The numerical simulations 
(Fig. 6) strongly support this result. 
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IV. SUMMARY 



In summary, we have considered slightly overloaded 
fiber bundles, and investigated how the fiber breaking 
rate progresses. It has a minimum after a number of 
steps to of the degradation process, and we have demon- 
strated that the total bundle collapse occurs near 2t - 
The demonstration has been performed for three differ- 
ent distributions of fiber thresholds, but the result is 
doubtlessly universal. Thus the fact that the breaking 
rate has a minimum predicts not only that a global fail- 
ure will occur, but also estimates when it will occur. 
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